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TECHNICAL NOTE HO. l8l6 


TRIANGULAR WINGS CAMBERED AND TWISTED TO SUPPORT 
SPECIFIED DISTRIBUTIONS OF L_LbT AT 


SUPERSONIC SPEEDS 
By Barrett S. Baldwin, Jr. 


SUMMARY 

The linearized supersonic flow theory has been applied to the 
calculation of the camber and twist required for various lift dis- 
tributions on a triangular wing with leading edges behind the Mach 
cone. 


Several specific examples have been derived and plots of the 
cambered surfaces, the corresponding lift distributions, and values 
of the theoretical drag coefficients are presented. The examples 
were chosen to avoid the singularities in lift distribution or 
camber that appear in previously given solutions and may be useful 
for practical applications since they avoid large pressure peaks 
and extreme angles of twist. 


INTRODUCTION 

The design qf cambered wings requires a knowledge of the rela- 
tion between the shape of the lifting surface and the distribution 
of lift over it. For airfoils of triangular plan form at supersonic 
speed, two examples of this relation are known; namely, the lift 
distribution for a flat wing at an angle of attack (reference l) and 
the camber required to support a uniform distribution of lift over 
the surface (reference 2) . 

The solution for the flat wing shows an infinite concentration 
of lift along the leading edges, while that for the uniform lift 
distribution shows that an extremely large angle of attack at the 
center line of the triangle would be required. Thus, both solutions 
show singularities which the designer may wish to avoid. It was 




2 


mca th Ho. 1816 


suggested to the present author by R. T. Jones that these undesir- 
able s ingular ities in the lift and camber distribution could be 
avoided, at a particular lift coefficient, by an appropriate choice 
of the functions used to represent the lift distributions. In the 
present report, several examples of such nonsingular lift distribu- 
tions and their corresponding camber shapes have been derived. 

Following a brief review of the method, a series of functions 
capable of representing a variety of lift distributions is presented. 
It is shown that the desired property of the lift distribution and 
camber shape can be achieved by an appropriate relation between the 
coefficients of this series, numerical examples showing lift 
distributions arid camber shapes which involve the first few terms of 
the general series are then given. 


METHOD 

t Conical Flow Theory 

As in the usual linearized theory (references 3 and. 4) , the 
boundary conditions are satisfied in the horizontal plane rather 
than in the act ual plane of the wing. The coor din ate system, as 
depicted in figure 1, is chosen so that the apex of the triangle is 
at the origin. The positive x axis extends back in the plane of 
symmetry of the triangle, the z direction is vertical, and the plane 
of the triangle is approximately horizontal. The symbols used are 
defined in Appendix A. 

Thf> fractional perturbation velocities are denoted by u, v, and 
w (i.e., fractions of stream velocity) in the x, y, and z directions, 
respectively. The difference in pressure coefficient between upper 
and lower surfaces is given by 



and the slope of the surface by 


dz 

dx 


= w 


The line arized equation is solved using the method described in 
reference 4. Solutions are found for a Mach number of */£, but may 
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"be extended to other Mach numbers by the use of the -well— known Prandtl 
transformation. 


In reference 4 it is shown that a special complex variable 
exists such that if f (£) is any analytic function of t . then 

u(£) = real part of f (£) and w (£) = real part of ^ — i 

is a conical solution of the linearized equation. 


The complex variable £ represents a transformation of the 
conical flow field which causes no distortion of the airfoil coordi- 
nates on the z=0 plane inside the Mach oone, but places the Mach 
cone along the real axis in the t, plane from £ = — <» to — 1 and 
from 1 to to. The exact expression for £ is , 


i = 


2e 

1+e 2 


where 


e _ y+iz 

x+ a/x 2 — y 2 — z 2 

The use of the variable £ results in a simplification of 
boundary conditions on the plane z=0, since it can be shown that 
£ approaches y + iz as z approaches zero, 
x 

To yield the lifting case, f(£) must be so chosen that its real 
part will vanish everywhere on the real axis in the £ plane except 
on the wing, which extends from £ = -m to £ = -hh. In addition, the 
real part of f(£) must be antisymmetric with respect to the real axis 
in the £ plane. If the wing lies entirely within the Mach cone, 
f ( £) must vanish at infinity. This latter requirement arises 
because in the £ plane the Mach cone is to be found at infinity in 
any direction as well as along the reel axis as previously described. 

For wings inside the Mach cone, the integration for w can be 
taken along the real axis from £* = 1 to £ * = £ making w zero 
everywhere on the Mach cone if f(£) is zero at infinity. The term 
£ * ip a dummy variable of integration replacing £ in the integrand. 

Thus, to represent the lifting surface and its pressure distri- 
bution, f (£) must have the following properties: 

1. f(£) must be an analytic function of £ . 
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2. f (£) must vanish at infinity. 

3. The real part of f(£) along the real axis must represent 
the desired lift distribution of the wing. 


4. The real part of f( t,) must he antisymmetric with respect to 
the real axis in the £ plane. 

When checking a prospective f(t, ) against these four require- 
ments, £ must he considered a complex variable. However, once 
these requirements are met, all subsequent calculations are made on 
the z=0 plane, which corresponds to the real axis in the £ plane, 
so that £ is equal to y/x, a real variable. Thus in the final 
expressions for u(£), w(£), and — (£) , £ is equal to y/x. 

An example of a function satisfying the foregoing requirement is 


Here 


f(£) 


m 


u(0 = for l?|Jm 

and 


u(£) = 0 for |£ J =? m 


or 


u(x,y,0) = 7 e=== over the wing. 

fZ 2 y 2 


Since 


w( 0 = real part 
it is found that 



w(£) = -jE (VI-® 2 ) for |C| =m 


and 
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*0 - f ° r Ul*» 

This solution gives the lift distribution for a flat triangle at an 
angle of attack. (See reference 1.) 

As a second example, choose 


f ( 0 = i log 

■where in this case it is found that 

u(£) = K for |£ | ^ m 

and 

u(£) =0 for |£ | ■t m 

Integration shows that 


v(£) = 


a/i-ei 2 

( cosh x 

l+m£ 

+ cosh 1 

1-mt; 

^ — cosh 1 

1 

m 


£+m 


/ m 

i 


for the whole region. This solution gives the required camber and 
twist for a uniform lift distribution. (See reference 2.) 


General Expressions for Lift, Drag, and Camber 
The camber is determined from the relations 

g- v( E ) -»(|) 

or 

< * ,7) ' j£i v (^) * 


(la) 
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The integration is taken from x* = |y| to x* = x so that the 
stream surface in "which the "wing lies may start at z=0 on the 
Mach cone. Thus, the shape of the surface is obtained in the form 
z = z(x,y). 

Since the wing surface is conical, its shape can be expressed 
in the form — = — (£), that is, — = function of £ where 

^2, Furthermore, the integration for z can be accomplished 

in terms of £, thus 



(lb) 


The magnitude of camber is related to the lift coefficient by 
finding the common factor of u, w, and — which will yield the 
desired lift coefficient. z 


m total lift . . 1 r TAP ^ 
q X wing area S Jg J 1 



u ds 


If u = Bu* where B is the desired factor. 


For the symmetrical triangular wing extending from x=0 to x=l 
along the center line and from y = — mtoy = + m along the trail- 
ing edge, the integration can be accomplished with a single integrar- 
tion in terms of £ . Since £ is constant in the element of area 
enclosed by a triangle the base of which is dy along the trailing 
edge and the vertex of which is at the vertex of the wing, 



-m 


and, for a symmetrical distribution. 
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CL = 



u*(0<iC or B = | X — |a— ^ 

J **(£)<*£ 


.( 2 ) 


The coefficient of drag 


due to lift is obtained as follows: 


C D “ 




u x wis 


The Integration over the wing can again he reduced to a single 
integration in terms of £ . 


f “(CMOn (3) 

Jo 


No additional term for leading— edge suction is required because 
the pressure distributions chosen have no infinite pressure peaks 
along the leading edges. 


SERIES DEVELOPMENT FOR AN ARBITRARY LHT DISTRIBUTION 
Elements Suitable for Series Representation 
The analytic functions 


fo( £) = i log ^ 
t — m 


w 


and 


tzn-l(0 = 


231-1 


n=l, 2 , 3 , etc. 


(5) 
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can "be superimposed to yield an arbitrary lift distribution extending 
from -m to +m. The function of fpn— i is similar to that given in 
reference 5* The exponent is restricted to odd integers so that the 
lift distribution will be an even function of £, or symmetric with 
respect to the center line of the triangle. Such a solution would be 
represented by a series of the form 

*«) - » i k g ♦ i 

n=i 


The velocity components corresponding to the individual terms are 


Usd-iCC) = real part 


Tzn-iU ) 


(6) 


and 


■» 2 n-i(£) = real part 



Vl -^ 2 


ar'an-i) 


( 7 ) 


It can be seen that this family of functions has the required 
properties of being analytic, vanishing at infinity, and having real 
parts which are antisymmetric with respect to the real axis in the 
£ plane. Also the uptwi terms are zero everywhere on the real axis 
except at | (; | = m. 

As an example of the requirement that the real part cf f (£) be 
antisymmetric with respect to the read axis, it will be shown that 
i (£— a/^ 2 — m 2 ) meets this requirement. For this purpose £ can be 
taken equal to since the argument can be restricted to the 

region where z approaches zero. It is Immediately apparent that 
the real part of i £ is anti sy mmetric with respect to the real axis. 
That the real part cf i ^ 2 -® 2 has this symmetry is not obvious but 
it can be seen as follows : 
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Let 


£+m = rie^^ 1 , £-m = 


•where ri is the length of the line Joining the point, -m, on the 
real axis to the point, £ , in the £ plane; and <pi is the angle 
this line makes with the positive direction along the real axis. 
Similarly, r 2 is the length of the line joining •a to { and <fe 
is the angle this line makes with the positive real axis. 

Then 


£ 2 -m 2 = (£+m) (£-m) 


= r x e 


Tze^ z 


ifcpi-Kpo) 

= rir 2 e 


i V t; 2 -® 2 = ±>JtiTz 


- i */ rxr 2 

and 

real part (i*/ ^-m 2 ) = — Vrir 2 sin 


From t he defi nitions of and <%, it follows that the real 
part of is always negative above the real axis and always 

positive below. This is seen by starting at a point above the real 
axis with cf>L and cfe less than jt and following a path which does 
not pass through that part of the real axis between +m and -m. 
Following a path that does pass between +m and -m is prohibited 
because the real part of f(£) will be discontinuous there if it is 
antisymmetric. The real part of f(£) is not discontinuous along the 
real axis outside of +m and -m because it passes through zero there . 
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Orthogonality Relations 

'The Upn — t terms are orthogonal in the interval -m to +m and 
are zero at the limit s of the interval. Any reasonable even conical 
lif t distribution can he expanded in an infinite series of Uq and 

Upr>— 1 . 

Let 


C = m sin 0 

for |£ | ^ m on the real axis, then 

fyn—i (m sin 0) = i(sin 0 — i cos e) 211 - 1 

= i 1_en+1 (cos 0 + i sin 0) 211-1 
- e i(2nr-i)0 

and finally 

fan-iCm sin 0) = (-l) 33 ^ 1 [cos(2n-l)0 + i sin(2ur-l)0 1(8) 

since 

Uyrv.1 (m sin 0) = real part [ far>-i (m sin 0)] 

uprwi fm sin 0) = ( — l) 21- 1 cos(2nr-l)0 (9) 


To establish the orthogonality of these functions, it is necessary 

n- 

to evaluate / 2 Uan-i d0: 


1 


7T ‘ TT 

2 Uptwi U 2JS-1 d 0 = 2 J z -008(2117-1)0 cos (2r— 1)0 d0 

~ 2 


or 


rZ 

I s 


2 — 1 upr- 1 30 — 


0, n f 


r u > r ^ 
t T> n = ri 


( 10 ) 


Equation (10) is the required orthogonality relation. 
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Let h(£) "be any even function which is zero at £ =± m. 
Expanding h(m sin 0) in the Fourier series. 


00 

h(m sin 0) = ^ h 2 nr-i (— l) 3 *"” 1 cos(2nr-l)0 


n=i 


and 


rb-1 ]|. 

b 2 n— i = (—1) — / h(m sin 0) cos(2n-l)0 d0 


( 11 ) 


The teim aouo = ao« can he added to h(m sin 0) to give a lift 
distribution which is not zero at £ = ± m. 

Then 


00 

u’(m sin 0) = a Q jr + ^ han-iC-l) 11-1 oos(2n-l)0 


( 12 ) 


n=i 


or 


00 

u '(£) = ao* + h2n— l Uart— i(£) 


( 13 ) 


n=i 


Since the real part of ^i log = rt for |£ |- m and the 


terms Uptwj are the real parts of 1 \ 


2 2 \ 2n “ 1 

1 it can he 


seen that the f*(£) which corresponds to u f (£) in equation ( 13 ) 
is 
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which, in terms of £ = m sin 0 is 


f*(m sin 0) = i a Q log 


sin 0+1 
sin 0—1 



n=i 


■^2n— l 


(-D 


H — 1 


[cos(2n-l)0 


+ i sin(2nr-l)0] 


(15) 


By the use of the orthogonality relation (equation (10)) it can 
he shown that no net lift is developed hy Use— i (£) when n is' 
greater than one. 

The total lift is proportional to 


J' J'u&s -f u(£)d£ 


then 


/m p 2 

Uptw- i ( £)d£ = / Ugn-i (m sin 0) m cos 0 


d0 


and 


n— 1 

= (—1) m / cos(2nr-l)0 cos 0 


d0 


Tt 




Uenr-i Hi 


nut, n=l 


So ^ (E,iE '{o;^i} 


(16) 
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Prom this relation the factor B, •which relates the magnitude of 
camber to the lift coefficient, can be expressed in terms of a Q 
and bi, the first two coefficients of the Fourier expansion of 

u*(0- 

Therefore, if 


then 


or 


00 

u (0 = B + ^ t>2nr-i 


n=i 


B = 


m 


Cl 


m 


CL 


**(£)&£ ^aoUW+bi^^ 


B = 


jt(4a 0 +bi) 


(IT) 


Condition for a Finite Talue of w(£) at £=0 
The expressions for Wan — 1 ani Wo all have the same type of 

_i -1 

infinity at £=0, namely, cosh By an appropriate relation 

between the coefficients of the Fourier series, the resultant 
coefficient of cosh -1 ^ may be made to vanish and the infinity 
in w( £) at £=0 will thereby be avoided. 


.According to A. Busemann in reference 3* the necessary and 
sufficient condition that there be no infinity in w at £=0 is 
that df=0 at £=0 . 

Differentiating equation (15) , 

(df) 0 =o = i £ 2ao+ Z (-l)^^-!) W-i J d 9 ( £=0 when 0=0 ) 
Therefore the necessary and sufficient condition that there be 
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no infinity in w at £=0 is 


2 & q + 




= 0 


n=i 


( 18 ) 


This relation restricts the lift distribution to the extent 
that the value of u at £=0 must be less than the average value, 
provided that no negative lift is allowed at any part of the wing. 
Also the value of u at £=m depends on the shape of the distri- 
bution of u over the whole wing. 


Ex a mples of Functions Contained in the Series Development 

Evaluation of w(£) and §■ (£) becomes increasingly complicated 
for higher values of n. The following is a list of the expressions 
for u(£)> w(^), and & (£) of the first three terms of the series: 
(They were derived by the use of equations (la), (lb), (^) , ( 5 )* (6), 
and (7). A sample derivation is given in Appendix B.) 


Uo(£) = * 


v 0 (C) 


¥( 6 ) 


UlU) 

vi(£) 


= a / 1 "”* 2 


^cosh - 1 


l+m£ 

£+m 


+ cosh' 


1 


^ £ cosh -i I 1 

I / m |£ 


a/l-^h 2 / £+m 


m m 


cosh' 


1 


l+m£ 

£+m 


- £ r^l) ‘ w 


_ Vm 2 -^ 


m 


= - */l -£ 2 - - cosh -1 I 7 1 + - [f fsin -1 ^ 

m m | SI * L \ / 

- e ( 3 ^- x * /E3 0 ] 
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3«> =5 lil +i 


K( Vl-m 2 ) — E( Vl-m 2 ) 


— — (oo) 

m ' 


Uq(C) = 


l^-m 2 


nT 




„ 3(?) - - 2 £ (X^) v= ♦ « » (sin- ^3S, ^3) 

— — F ( sin- 1 - ^S , */l^F V 4 C V?W + 2 coatr 1 ! 

^(0 = -|(a/hF- Ul cos-h) -£[ (2H 2 ) - 3 |£|cos-h] 

+ ^2- E ( Vl-m 2 ) - - K ( ■*/l-ffl 2 ) + - (o») 
m mm 


The preceding egression for ^ (£) is derived from an integra- 
tion of only those terms of w which, are variable over the range 
Is I ^ m plus a constant term which is needed along the line, £=0. 

For each there will be one term from the integration of 

cosh 1 jk, which iff infinite at £= 0. When the individual solutions 
b — "| 

are superimposed in a series, the cosh =• terms are removed by 

imposing equation (l8). Terms of ^ which are to be removed by 

equation (l8) have been denoted by the symbol («°) in the table. 


The elliptic in tegrals F and E are complete and hence constant 
for |? | ^m. Also */£2-m 2 is taken to be zero for |?|- m, since w 
was to be only the real part of the function resulting from an inte- 
gration. 

An arbitrary dihedral angle can be added to the wing surface 
along the wing center line without changing the lift distribution as 
long as no part of the wing is far from the z=0 plane. Thus a term, 
a constant timffs t, can be added to the equation of the surface 
without changing w (£) oru((). The value of this constant determines 
the dihedral angle along the wing center line. It can be seen 


H | 
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mathematically that a term £ equals constant times £ will he 
proportional to and w (which equals “) will he zero. 

Such a term added to the equation of the wing surface has the 
effect of shifting the stream sheet in which the wing lies so 
that it does not pass through z=0 on the Mach cone. However, 
such a departure iB allowable under the condition that no part of 
the wing surface sheet lies far from the z=0 plane . This 
represents an approximation which is valid to the same order of 

magnitude as the practice of calculating ^ in the z=0 plane 

rather than in the wing surface. For example, Stewart* s solution 
for the flat triangle yields a V-shaped wing-surface cross section 
which is equivalent to a flat triangle when the dihedral angle is 
removed. 


Examples of Specific Lift Distributions 

Four cases have been worked out using f 0 (£), f (£), and f (£) 
in various combinations which comply with equation (18) . The 

expressions for f (£), u(£), w(£), and ^(£) are listed in each case. 

Graphs are drawn of u and ^ for Cj=1.0, nt=0.577 ia figure 2. 

For purposes of comparison, the flat triangle and constant load 
triangle are also included in the figure. 

The wing-surface shapes of figure 2 are drawn to a square scale. 
(The unit of length in the horizontal direction is equal to the 
vertical unit of length.) The z=0 plane is the horizontal line 
above the wing. The semi span of the wing is in each case equal to 
0.577 of the radius of the Mach circle. The leading point of the 
wing would be at the origin of the coordinate system, a distance 
ahead of the plane depicted in the figure equal to the. radius of the 
Mach circle in that plane. 

The cases considered are most easily evaluated by the use of 
the following tabulated expressions for u, w, and 


00 

u(£) = B |a 0 u 0 i-^bai-a. uair-i (£) 
npjL 
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v(t) - B j^a 0 Vo + ^2Qr-l V 2 lt-1 
n=i 

|U) = B [ao— + J Wa (g)J 

D=1 

The term bp^ will "be zero for n>e in the oases considered, 
here. 

The coefficient a 0 may- he taken equal to- 1 or 0, 3f a Q is 
equal to 0, h may he taken as 1, and. h s can he chosen with a 
value such that the cosh - * 1 i term does not appear in w T ( t) or the 

(op) term does not appear in ^ (£), This choice is equivalent to 

JU 

satisfying equation (l8). 

From equation (17), 

* °±— 

*( 1 * ft o +lx i) 

/ 

The coefficient of drag due to lift was found hy a graphical 
integration of equation (3) # 

A term D|£ | is added to £(£) to make allowance for the 

arbitrary dihedral angle at the center line of the triangle. For 
example, D can he given the value required to place the leading 
edges at the same vertical height as the center line. This value 
was used in the cases considered here. 
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Case I 


" 1 1 [ lo e IS - 1 

u(o =i (*-i 




cosh. 1 I I + cosh" 
I t+m I 


I •/ m 


- 1 Tk ( - E ( ./E5P) | 


i(a=aj (s& COST' |i! 2 L| - £= cosh"* |£=£|) 

x' =>' 2it \ m \ m ) £ +m “• I C - ® 1 ' ' 


— £ ^*/l— £ 2 — |C| 008 1 ^) ” § [k ( ^l - ™ 2 ) — E ( Vi-m 2 ) J “ D |d ^ 


Cl - ! 


m = 0.577 


z 

X 


D = 1.486 

£ 


CD = 0.25 


u 


.05 
.10 
.15 
.20 
.22 
.24 
.26 
.28 
.30 
• 35 
.40 
.45 
.50 
.55 
.577 


-0.459 

-.429 

-.4o4 

-.383 

-.366 

-.361 

-.356 

-.353 

-.350 

-.347 

-.346 

-.350 

-.362 

-.384 

-.421 

-.459 


0 


.1 

.2 

.3 

.4 

.5 

.577 


0.182 

.188 

.201 

.228 

.271 

.341 

.500 
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Case II 



£+m 2 

l0g f^ + 3^ 


(£ 


«« ■ e 


« + jTf (4£ 2 -“ 2 ) y» 2 -£ ! 
om 


*<0 - 


2 r 


cosh 


l+mg 

£-Hm 


+ cosh 1 


l-m£ 

t,-m 



_ _§_ + -2_ 
3m s 3m s 


( 4 + m 2 ) E ( //l-m 2 ) — 5m 2 E 



§(» = 

U--C 2 - Ul oos-*^ 


1-Hn£ 


£+m 


C-m 

m 


cosh' 


l 


l-m£ 


) 


2 

m 


4 

a? 


( 2 +^y^-P - 3 U 1 003-1 $ 


+ 



(4+m 2 )E( ,/l-m 2 ) — 5m 2 £(,/ 1-m 2 ) 



C 


L 


1 


m = 0.577 


D = 22.481 


Cjj = 0.26 


X_ 

.05 

.10 

.15 

.20 

.22 

.24 

.26 

.28 

.30 

.35 

.40 

.45 

.50 

•55 

.577 


z. 

X 

i 

M 

-0.546 

0 

0.196 

-.489 

.1 

.204 

-.442 

.2 

.224 

— .406 

.3 

.253 

-.380 

.4 

.286 

-.373 

•5 

.302 

-367 

-.363 

-.361 

-.360 

-.365 

-.382 

-.410 

-.449 

-.504 

-.546 

.577 

.250 
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Case IH 

f(£) = ^ *[; <e-«/F^) * ~ 


' *?{ S 3 (1 -5 2,s/a + ^[(2-a 2 ) e (7E5?) e (^5?)]! 


£ 

2 


(£) - F 1 f ^» [<*•£*> A? - 3 IE | oob-V 

+ E (v'l-n?) -II s E (Vl-sn 2 ) -B|J| j. 


°I- 1 


E 

0 

.05 

.10 

.15 

.20 

.22 

.24 

.26 

.28 

.30 

.35 
.4o 
.45 
.50 
.55 
• 577 


m = 0.577 


z 

X 


-0.633 

-.549 

-.481 

-.429 

-.394 

-.385 

-.378 

-.373 

-.372 

-.372 

-.385 

-.414 

—•^57 

-.515 

-.580 

-.633 


D = 10.50 


i 

u 

0 

0.212 

.1 

.224 

.2 

.242 

.3 

.279 

.4 

.300 

.5 

.265 

.577 

0 


Cjp 0.28 
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Case IV 


f(?) - 1 1 [ lo8 f!= ■ - 1 < « + 3? « -J^y 


n(t) = Si 
3* 


3m 


w(C) 


= °L f 


3* 


cosh. 


—l 


1+mC 


£+m 


+ cosh 


—l 


i-*E| 






+ rra [(l+m 2 ) I (VT 5 P) - 2 m 2 E (/EiF) ] l 

3^ L J J 


z ( ¥ \ J a/i-di 2 /d®. 

? l d = 3«t — 5”V~ 


cosh 


—l 


l+m£ £-m —i 

715 5“ Gosl1 


l-m£ 

£-5 


) 


- - ( J'L-i* - Id cos' 1 0 [(2+t 2 )/K 2 - 3 Id COS 

m 3m 1 1 


£] 


+ ^3 [(W) E (7 1-m 2 ) -^m 2 ^: (/ 1-m 2 )] -D |d } 


°L " 1 


m = 0.577 


D = 11.984 Cj) = 0.26 


Z 


JL_ 

I 

i 

i 

-0.517 

0.28 

.05 

-.469 

• 30 

.10 

-.430 

• 35 

.15 

-.398 

.40 

.20 

-.376 

.45 

.22 

-.369 

• 50 

.24 

-.364 

• 55 

.26 

-.360 

• 577 


z 


I 

t 

U 

-0.357 

0 

0.192 

-.356 

.1 

.198 

-.359 

.2 

.216 

-.371 

-3 

.245 

-.39^ 

.4 

.281 

-.428 

• 5 

.316 

-.476 

-.517 

• 577 

• 333 
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EVALUATION OF RESULTS 

As can "be seen in figure 2, the theoretical coefficients of 
drag due to lift of the camber shapes considered here are all larger 
than that of the flat triangle. In practice, however, it is possible 
that case II and case IV would have lower drag than the flat triangle 
at some values of lift coefficient because of the suppression of 
large adverse pressure gradients. 

Quite possibly case II and case TV would be capable of attain- 
ing a higher lift coefficient than the flat triangle because of 
their favorable pressure gradients, although the flat triangle 
pressure distribution would be superimposed' on their pressures at 
lift coefficients other than that for which they are designed. 

Case HI has the most radical camber, which might cause 
uncalculated additions to the drag as well as flow separation. 

Theoretically the dihedral angle in the surface shape is 
arbitrary. In practice, however, variation of this parameter might 
result in an optimum of characteristics at other angles than those 
chosen in figure 2. 

While the linear theory cannot be expected to apply for values 
of lift coefficient greater than about 0.25, the calculations were 
made at a lift coefficient of one to better Illustrate the wing- 
surface shapes. This value of lift coefficient, 0.25, is also an 
approximate limit beyond which the assumption that no part of the 
wing lies far from the z = 0 plane is invalid. 

These solutions are applicable to swept-back or other wings 
obtained by cutting out the rear of the triangle along lines such 
that no part of the remaining wing lies within the zone of influence 
of the parts removed. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX A 
Symbols 

common factor of u, w, and z which relates the magnitude of 
camber to the lift coefficient 


B 
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°D 

C L 

D 

E 

€ 

fig) 

E 

K 

m 

M 

<1 

<1 

3 

u,v,v 


coefficient of drag due to lift 
lift coefficient 

constant, the value of which determines the arbitrary 
dihedral angle at the center line of the triangle 

elliptic integral of the second kind (either complete of 
incomplete) 

y+iz 

complex variable equal to — : — ■■■■ ■■ — 

x+ v^ 2 -y 2 -z 2 

analytic function of the complex variable £ such that the 
real part of f(£) equals u 

incomplete elliptic integral of the first kind ' 

complete elliptic integral of the first kind 

horizontal slope of the leading edge of the wing 

Mach number (equal to J2. ) 

the difference in pressure coefficient between upper and 
lower wing surface 



wing area 

fractional perturbation velocities as fractions of the stream 
velocity 


u* u divided by B 

V stream velocity 


x,y,z Cartesian coordinates 


g 


equal to 


also used as complex variable equal to 


2e 

1+e 2 


£ f ,x f dummy variables used In integrations when a limit of integra- 
tion is t or x 
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absolute value 

terms in the equation of the ving surface vhich are to he 
removed by the superposition of other independent solutions 


APPENDIX B 

SAMPLE CALCULATION OF EQUATIONS IN TABLE OF u, v, AND — 


From equation (5) 




Since the subsequent calculations are made in the plane 
= 0, £ is a real variable equal to %■. 


df i = 


= i(^i \i 


From equation (7) 


r 1 

vx = real part — i — df f x 

Vl = i 

iJl t' “A /£‘2 h»? 

a S' * "A- E 2 - oosh -1 i 

/’‘ffi- df = d£ ♦ l* ^ df 

^ V £ f 2 -m? « ^ C 2 -® 2 o m , 2 -m? 


d p.r&L d£ for £<m 


The last equality holds because the second term on the right 


The last equality holds because th 
of equation (B3j is Imaginary for £< m 
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With the substitution, r 2 = 


Jdtf 

1-ffl 2 



r g dr 

//L-r 2 a/L— ( l-m2)r 2 



l-g]2) r g dr 

,/l-r 2 



dr 


1-r 2 a/L— ( 1-ni 2 ) r 2 


The two terms on the right of the last equality are recognized as 
complete elliptic integrals of the second and first kind# respectively, 
so that 


J x Jv^ 


d£ ? = EC^/L-m 2 ) - for 




(b4) 


The same substitution leads to the result 



These expressions all apply for £ negative as well as 
positive because wi is an even function about £ =0. 

Thus from (Bl), (B2), (B3), (B4), and (B5) 


w x 


= £ v^— C 2 — cosh 1 — — E C/l-m 2 ) + J- for | t ) < m 
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and 


Vi 


- i W 2 - 


cosh 


l 


E 


(■ 


sin” 1 

Vl-m 2 y 


+ F 


(sin -1 \ 1 for |£|>m 

V Vl-m2 y J 


From equation (lb) 


Zi 

x 




icl 


viCC*) 

c 2 


dC 


(B6) 


Since the dihedral angle at the center of the triangle is 
arbitrary, the elliptic integrals, vhich are complete and hence 
constant for |£[ = m, need not be included in the integration. 
However, to achieve the correct angle of attack at the center line, 
they must be taken into account. This operation can be accomplished 
most simply by the use of equation (la), which is 


‘fw v (£) 


dx f 


Thus along the center line, the term of 
integrals in the expression for Vi is 


zi from the elliptic 


P X - K( Vl5) - E( Vl-® 2 ) 
J a m L 


dx» = J 


K — E 


(B7) 


The terms resulting from the integration of cosh -1 |y| will 
he removed by a superposition of other solutions according to 
equation (l8) to avoid an infinite angle of attack at the center 
line. Such terms are denoted by the symbol (“>) in the expression 

for 


The only term of Vi not yet considered is ^ J 1— t, z . Inserting 
this term in (b 6) yields 

1^1 1 V-i-i - at, = 1 

m fcp ® m 




i 

t a 


yiT 5 - it i cos- 1 1 


(B8) 
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Thus from (b 6), (B7), and (B8) 


Zl 

X 


1 

M 


- Ul cos” 1 5 



m 


K( vGF ) - E( ) 
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Figure /.- Cambered triangular wing at an angle 
of ottack showing the coordinate system. 
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Figure 2.- Cross sections of wing surface shapes and lift 
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distributions (Ci_=i). 
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